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ABSTRACT 

In evolved supernova remnants (SNRs) interacting with molecular clouds, a highly inhomogeneous structure 
consisting of a forward shock of moderate Mach number, a cooling layer, a dense radiative shell and an inte- 
rior region filled with hot tenuous plasma is expected. We present a model of nonthermal electron injection, 
acceleration and propagation in that environment and find that these SNRs are efficient electron accelerators and 
sources of hard X- and 7-ray emission. A forward shock of velocity v, > 100 km s" 1 with an ionized precursor 
propagating into the molecular cloud accompanied by magnetohydrodynamic turbulence provides a spatially in- 
homogeneous distribution of nonthermal electrons. The energy spectrum of the nonthermal electrons is shaped 
by the joint action of first and second order Fermi acceleration in a turbulent plasma with substantial Coulomb 
losses. Bremsstrahlung, synchrotron, and inverse Compton radiation of the nonthermal electrons produce mul- 
tiwavelength photon spectra in quantitative agreement with the radio and the hard emission observed by ASCA 
and EGRET from IC 443. We distinguish interclump shock wave emission from molecular clump shock wave 
emission; particles reach higher energies in the interclump shock and that is the likely source of 7-ray emission 
and radio synchrotron emission. Spatially resolved X- and 7-ray spectra from the supernova remnants IC 443, 
W44, and 3C391 as might be observed with BeppoSAX, Chandra XRO, XMM, INTEGRAL and GLAST would 
distinguish the contribution of the energetic lepton component to the 7-rays observed by EGRET, constraining the 
cosmic ray nuclear component spectra in these SNRs. These data would provide a valuable tool for studying the 
complex structure of molecular clouds where SNR radiative shocks interact with dense molecular clumps. 

Subject headings: acceleration of particles — cosmic rays — ISM: Supernova Remnants — radiation 
mechanisms: nonthermal 



1. INTRODUCTION 

Supernovae and their remnants have long been identified as 
likely sites of galactic cosmic ray acceleration. Evidence for 
electron acceleration first came from observations of radio syn- 
chrotron radiation. The observation of nonthermal X-ray emis- 
sion from SN 1006 by ASCA has provided evidence for electron 
acceleration up to ~ 100 TeV (Koyama et al. 1995; Reynolds 
1998). The proton cosmic ray component is more difficult to 
observe, but it has long been recognized that pion decays from 
collisions with interstellar gas could give an observable flux 
of 0.1-1 GeV photons. Initial estimates of 7-ray emission 
from supernova remnants (SNRs), concentrating on the com- 
ponent resulting from pion decays, were made by Chevalier 
(1977), Blandford & Cowie (1982), and Drury, Aharonian, & 
Volk (1994). When 7-ray emission was apparently detected 
from SNRs by CGRO, the 7-ray spectrum could not be fitted 
by a pure pion decay spectrum and some other component was 
needed (Esposito et al. 1996; Stumer et al. 1997; Gaisser, 
Protheroe, & Stanev 1998). In addition to the pion decays, the 
relevant processes are bremsstrahlung emission of relativistic 
electrons and inverse Compton emission. Gaisser et al. (1998) 
modeled these processes in detail in order to fit the observed 
7-ray spectra of the remnants IC 443 and 7 Cygni. They as- 
sumed acceleration to a power law spectrum in the shock front 
and determined the spectral index, electron to proton ratio, and 
the upper energy cutoff. 

The evolution of the relativistic component in a supernova 
remnant has been modeled by a number of groups. Sturner et 
al. (1997) assumed particle acceleration in the shock front with 



an E~ 2 energy spectrum and followed the evolution of the par- 
ticle spectrum. Their SNR model was homogeneous: a factor 
4 density jump at the shock front, constant density in the inte- 
rior, and a constant magnetic field in the interior. Their shock 
acceleration model had a diffusion length parameter that was 
important for the maximum energy reached by the particles. 
They modeled the nonthermal emission from IC 443 as syn- 
chrotron emission in the radio and bremsstrahlung in 7-rays. 
de Jager & Mastichiadis (1997) dealt with the same processes 
in W44, as well as inverse Compton emission. They noted that 
the observed radio spectrum is flatter than would be expected 
from shock acceleration of newly injected particles and sug- 
gested that the particles originated from a pulsar in the super- 
nova remnant. Recently, Ostrowski (1999) showed that second- 
order electron acceleration by the turbulent medium just after 
the shock could flatten the spectra to account for the radio ob- 
servations of IC 443. Baring et al. (1999) presented calcula- 
tions of the broad-band emission from nonlinear shock models 
of shell-type SNRs. They used Sedov adiabatic shock dynam- 
ics in a homogeneous medium and a Monte Carlo simulation of 
the particle acceleration taking into account the nonlinear shock 
structure. The set of models considered by Baring et al. (1999) 
covers the range of shock speeds 490 < V5 < 4000 km s" 1 and 
ambient medium number densities 10~ 3 < n < 1 cm" 3 . 

Massive stars that are the likely progenitors of core collapse 
supernovae are expected to be spatially correlated with molec- 
ular clouds. The remnants that are likely to be 7-ray sources 
in CGRO observations (Esposito et al. 1996) also show evi- 
dence for interaction with molecular gas. Chevalier (1999) re- 
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cently studied the evolution of supernova remnants in molecular 
clouds and concluded that many aspects of the multiwavelength 
observations could be understood in a model where the rem- 
nants evolve in the interclump medium of a molecular cloud, 
which has density of 5-25 H atoms cm" 3 , and become radia- 
tive at radii <~ 6 pc. The compression in the radiative shell is 
limited by the magnetic field. Molecular emission occurs when 
the radiative shell collides with molecular clumps. IC 443 is 
the remnant with the best evidence for high energy emission 
and also shows clear evidence for interaction of the shock with 
a molecular cloud (e.g., Burton et al. 1990; van Dishoeck et 
al. 1993; Cesarsky et al. 1999). 

Most of the current models of particle acceleration in SNRs 
have dealt with adiabatic remnants in a homogeneous medium 
(see, however, the model of Boulares & Cox 1988 for the 
Cygnus Loop and the model of Jones & Kang 1993). In con- 
trast, our aims here are to examine in detail the nonthermal 
emission of radiative supernova remnants in molecular clouds 
using a kinetic model of electron injection, acceleration and 
propagation, including the inhomogeneous structure deduced 
from multiwavelength SNR observations. We use global MHD 
(magnetohydrodynamic) models of supernova remnants inter- 
acting with molecular clouds as described by Chevalier (1999), 
including the theory of radiative shock structure by Shull & Mc- 
Kee (1979) and the kinetic model of electron injection and ac- 
celeration by shocks from Bykov & Uvarov (1999). The model 
described here is relevant to the mixed-morphology SNRs (Rho 
& Petre 1998), if they are interacting with molecular clouds. 
These remnants comprise a sizable fraction of the Galactic SNR 
population. 

The plan of our paper is as follows. In § 2, we treat the en- 
ergy spectrum of energetic electrons as they pass through the 
regions associated with a radiative shock wave. We include the 
shock acceleration of particles from the thermal pool. In § 3, 
we present calculations of the nonthermal emission from the 
relativistic particles, including emission from shocks in clumps 
as well as emission from the radiative shock wave in the inter- 
clump region. Our model is specifically applied to the well- 
observed remnant IC 443, although we expect it to more gen- 
erally apply to other supernova remnants in molecular clouds. 
We also examine the ionization by energetic particles and the 
energy involved in the relativistic particles. A discussion of re- 
sults and future propects is in § 4. A detailed discussion of ex- 
pected range of densities, shock velocities and magnetic fields 
determined from multiwavelength observations of IC 443, W44 
and 3C391 was given in the paper by Chevalier (1999, and ref- 
erences therein). 

2. EVOLUTION OF NONTHERMAL ELECTRONS 

The structure of the flow expected in a radiative shock can be 
divided into the following regions (see, e.g., Draine & McKee 
1993): 

(7) preshock gas ionized and heated by UV (ultraviolet) emis- 
sion and fast particles from the shock; 
(77) shock transition region; 
(777) postshock cooling layer; 

(IV) cold shell of swept-up gas; 

(V) hot low density gas interior to the shock. 

Radiative shocks are subject to thermal instabilities (Chevalier 
& Imamura 1982; Bertschinger 1986) as well as dynamical 
instabilities (Vishniac 1983; Blondin et al. 1998). Although 
the effect of inhomogeneous magnetic fields on the stability of 
radiative shock flow is yet to be studied, MHD turbulence is 



expected in such a system. The transport and acceleration of 
nonthermal particles in the violent SNR environment are gov- 
erned by the MHD flow and depend on the ionization structure 
and on the spectrum of MHD turbulence. Shocks with veloc- 
ity V57 > 1.1 generate sufficient ultraviolet UV radiation for full 
pre-ionization of hydrogen and helium (He + ) in the preshock 
region, where the gas temperature is typically ~ 10 4 K (Shull 
& McKee 1979). Here, v S7 is the shock velocity measured in 
units of 100 km s" 1 , and n\ is the ion number density measured 
in units of 10 cm -3 . Shocks with v$i ~ 1.5 propagating in a 
medium with n\ <~ 1 photoionize gas ahead of the viscous jump 
to a depth TV, > 3 x 10 17 cm -2 . The shock transition in such a 
case is collisionless and supercritical unless the preshock mag- 
netic field is above — 4x 10" 5 n" 5 G. 

The geometry of the magnetic field is important for the ef- 
ficiency of electron injection and acceleration. It has been ar- 
gued (e.g., Levinson 1996; McClements et al. 1997) that quasi- 
perpendicular shocks are an efficient source of freshly injected 

electrons if their speed v, > (m e /m p )°' 5 c; the injection effi- 
ciency drops sharply below that threshold. On the other hand, 
Bykov & Uvarov (1999) have shown that quasi-parallel MHD 

shocks of moderate Alfven Mach number M a < (m e /m p ) °' 5 
are suitable sites of electron injection and acceleration. The 
SNR shock waves in molecular clouds are expected to be in 
this regime of Mach number. Thus we shall discuss below the 
quasi-parallel portion of the shock (magnetic field inclinations 
< 7r/4) propagating into a molecular cloud with v S i ~ 1.5. 

Nonresonant interactions of the electrons with fluctuations, 
generated by kinetic instabilities of the ions in the transition re- 
gion of a quasi-parallel supercritical shock, play the main role 
in heating and preacceleration of the electrons. The electron 
heating and preacceleration (injection) occur in the collision- 
less shock front transition region on a scale A of several hun- 
dred inertial lengths of the ions, (Bykov & Uvarov 1999). 
Here U = c/uj p i ~9x 10 6 «7°' 5 cm an d <-<Jpi is the ion plasma fre- 
quency. For a SNR shock in a molecular cloud, the scale A ~ 
10 9 cm is much shorter then the scales of all the other regions. 

To accelerate electrons injected in the shock transition re- 
gion to relativistic energy E(p), MHD turbulence should fill 
the acceleration region, which has a scale l a <~ k m (p) /v s . Here, 
k m = max-jTc,}, where ki(p) is the diffusion coefficient of an elec- 
tron of momentum p and ;' =(I), (III) refers to the regions defined 
above. To provide the conditions for efficient transformation of 
the MHD flow power to the accelerated electrons, the MHD tur- 
bulent fluctuation spectrum in the vicinity of the shock should 
extend to scales I > r e (p) (resonant scattering), where r e (p) is 
the gyroradius of an electron with momentum p. The MHD 
fluctuations responsible for the electron scattering are collision- 
less for particles with energies up to a GeV because the ion- 
neutral collision length is above 10 14 cm for typical parameters 
of the radiative shock described below (a ke V electron has r e w 
10 8 [B/10 _6 Gr 1 cm). 

To calculate the spectra of nonthermal electrons in the re- 
gions i = I - IV, we use a kinetic equation for the nearly- 
isotropic distribution function Ni(z,p,t): 

d d v d f d \ 

ki{p) + J_ * p 2 D i (p) f* + ^ * [p 2 Li(p)Nftl) 

oz l p l op op p L op 

The Fokker-Planck type equation (1) takes into account diffu- 
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sion and advection (with bulk velocity m,[z]) of an electron in 
phase space due to particle interaction with MHD waves and the 
large scale MHD flow (Bykov & Toptygin 1993). Here, U(p) 
is the momentum loss rate of an electron due to Coulomb col- 
lisions in a partially ionized plasma (e.g., Ginzburg 1979). The 
momentum diffusion coefficient D{p) is responsible for second 
order Fermi acceleration, and kj(p) is the fast particle spatial 
diffusion coefficient. For low energy electrons (i.e., E < Ec), 
Coulomb and ionization losses are important in the regions (7 
- V), except for the narrow shock transition region (77) where 
acceleration is fast enough to overcome losses and nonthermal 
electron injection occurs. The characteristic energy E c depends 
on the plasma density and ionization state, magnetic field, and 
fast particle diffusion coefficients h(p) and Dj(p). The corre- 
sponding momentum /?„ can be estimated as the point where 
the last two terms in equation (1) are equal. We used the en- 
ergy E c calculated simultaneously with the electron distribu- 
tion function as a convenient parameter to distinguish between 
different possible models of MHD turbulence in the postshock 
cooling layer of the radiative shock structure (see below in this 
section). 

The diffusion coefficients ki(p) and Dj(p) depend on the 
spectrum of the collisionless MHD turbulence, which is poorly 
understood. A plausible approximation for the diffusion coef- 
ficients ki(p) was assumed here. For all of the regions we used 
the following parameterization: 

{1, Pt<P<P* 
vp a /v*pl, p*<p<p** (2) 

cp 2 /(pl~fpy*) P**<P< Pm- 

Here, pj is the momentum of the upstream thermal electrons, 
p m is the momentum corresponding to the upper cut-off energy 
E m (see equation [3]), and p* and p** are defined in the next 
paragraph. The standard relation for the momentum diffusion 
coefficient A(p) = p 2 w 2 /9kj(p) (e.g., Berezinsky et al. 1990) 
was used. In the low energy regime p < p*, the particle trans- 
port is dominated by large scale turbulent advection (Bykov & 
Toptygin 1993). The large scale turbulence is due to stochastic 
MHD plasma motions on scales A longer than the particle mean 
free path due to resonant scatterings. The sources of MHD 
turbulent motions are the shock wave instabilities mentioned 
at the beginning of this section. The large scale vortex rms- 
velocity w, typically is a fraction of bulk speed m,. The spatial 
diffusion is energy independent below p* and kg) ~ w,-A. For 
p* < p < p**, electrons are scattered by resonant MHD waves 
presumably generated by the streaming instability of shock ac- 
celerated particles (e.g., Blandford & Eichler 1987; Jones & El- 
lison 1991). In this case w, is close to the Alfven velocity v a and 
ki(p) is the electron diffusion coefficient due to resonant wave- 
particle interactions. If the power spectrum of magnetic field 
fluctuations of Alfvenic turbulence for the resonant fluctuations 
is approximated as dBf/dk <~ k~ e then the power-law index a 
(equation [2]) is related to the index 8 as a = 2-9. An im- 
portant particular case corresponding to a = 1 can be written as 
kiip) = (l/3)r]ivr e (p). The parameter 77/ determines the scatter- 
ing "strength" and strong scattering (Bohm limit) corresponds 
to rji <~ 1. The diffusion model with a = 1 was successfully 
used for modeling the observed anomalous cosmic ray proton 
fluxes in the interplanetary medium (Ellison et al. 1999) and 
in the Monte Carlo simulations of nonlinear shocks in the shell 
type SNRs (e.g., Baring et al. 1999). The momentum p* in our 
model can be estimated from the equation (1 /3)rjivr e (p*) = w,A. 



MHD waves in the partially ionized plasma of the radiative 
shock flow are subject to ion-neutral damping (e.g., Volk et 
al. 1981). We define the momentum from the condi- 
tion r e (p**) = A,„, where A,„ is determined from the relation 
A,„ = 2wv a /Ti n and r,„ is the Alfven wave damping rate due to 
ion-neutral collisions (Kulsrud & Cesarsky 1971). For p > 
electrons are scattered by nonresonant small scale waves. 

The limits imposed on the diffusive shock acceleration of 
particles due to ion-neutral MHD wave damping have been dis- 
cussed by Draine & McKee (1993) and Drury et al. (1996). 
They considered the spectra of Alfven waves driven by the 
instabilities of shock accelerated ions in a partially ionized 
plasma. Drury et al. (1996) obtained the upper cut-off energy 
of a proton due to ion-neutral wave damping 

E m < v^ 7 r 4 - 04 n- 1 nf 5 GeV, (3) 

where n„ is the neutral particle density, n, is the ion density 
(both are measured in cm" 3 ), and VE\ j s tne tota j particle pres- 
sure normalized to 10% of the shock ram pressure. The ioniza- 
tion structure of the preshock region is inhomogeneous so the 
ion number density, temperature and other numbers in equation 
(3) must be taken at an upstream distance <~ ki(E m ) /v s from the 
collisionless shock transition. Application to our model yields 
E m ~ 10 GeV if vy? ~ (1.1-1.5) and ki(E m ) ~ 10 24 cm 2 s" 1 be- 
cause of the high ionization in the preshock region (I) and the 
postshock cooling layer (III) provided by the UV radiation. One 
can justify the same value of E m for the electrons because the 
synchrotron and inverse Compton losses are relatively unim- 
portant for E m ~ 10 GeV. 

Region I. In the preshock region of the shock with v S7 > 1.1, 
we assume a highly ionized plasma of temperature T4 ~ 1 up 
to depths ~ a few times 10 17 cm -2 upstream. For the IC 443 
preshock region, we adopted a density <~ 25 cm" 3 and v$i ~ 
1.5 (Fesen & Kirshner 1980). The magnetic field fly w 10" 5 G 
(Chevalier 1999), while B±^5x 10" 6 G (the inclination is as- 
sumed to be ~ 7r/6). The Alfven velocity here is about 5 km 
s" 1 . We assume a conservative diffusion model in the shock up- 
stream region where kj(p) has a = 1 and the moderate scattering 
"strength" 10 < 77/ < 100 (Ellison et al. 1999). 

Region II. In the shock transition region, injection and heat- 
ing of the electrons occur due to nonresonant interactions with 
strong MHD fluctuations generated by the ions. Following the 
model developed by Bykov & Uvarov (1999) we take 

**»»'*{%)'<%)'(;%)• < 4 > 

where C <~ 1 and SB > Bo. The Alfvenic Mach number of the 
shock, M a > 10, is moderate for our situation. Nevertheless, 
the effect of the accelerated particle pressure on the shock struc- 
ture could be important (see, e.g., Jones & Ellison 1991; Baring 
et al. 1999). An exact treatment of the effect is not feasible in 
our model because it would require a kinetic description of the 
ion injection. We accounted for the effect by correction of the 
total shock compression ratio S, and by introducing an extended 
modified shock structure with a precursor and a subshock with 
compression ratio S s < S t (e.g., Jones & Ellison 1991). In our 
case with Mach numbers M > 15, the estimated extra com- 
pression was about 10% and thus S, s=a 4.4. A subshock com- 
pression ratio (5. s m 3.3 is consistent with the assumption of sub- 
stantial heating of the precursor gas due to wave dissipation 
(Berezhko & Ellison 1999). These numbers are in agreement 
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with those from a hybrid simulation of nonlinear shock struc- 
ture (Giacalone et al. 1997). We found that the resulting high 
energy emission spectra are sensitive to the particular choice of 
precursor flow structure in the MeV regime. The keV and GeV 
emission is not sensitive to the structure if the compression ra- 
tios 5 t and 6 S are fixed. 

Region III. The postshock cooling layer has about a column 
density ~ 3 x 10 17 cm" 2 of highly ionized plasma with an initial 
density about S, times that of the preshock one. The Alfven ve- 
locity in the highly ionized portion of the region (III) is similar 
to that at the preshock region, <~ 5 km s" 1 . Large scale vorticity 
with an amplitude < 20 km s -1 on scales A < 10 14 cm may 
occur here because of shell instabilities of the radiative shock 
with v S7 ~ 1.5 (see, e.g., section 5 in the review by Draine & 
McKee 1993). This would dominate the low energy electron 
propagation, providing a spatial diffusion coefficient of <~ 10 20 
cm 2 s" 1 for the radiative shock parameters described above. 

Coulomb losses are important here for electrons with ener- 
gies below E c ~ 20 keV if the large scale turbulent velocity has 
a substantial longitudinal component with wni ~ 20 km s" 1 . 
Ec ^ 8 MeV is expected in the absence of a longitudinal com- 
ponent of large scale turbulence in the postshock cooling layer 
because the Coulomb losses are overcome by resonant interac- 
tion with Alfven waves. Since the Alfven velocity is relatively 
low (<~ 5 km s" 1 ), the value of Eq is much higher in that case. 
An accurate description of MHD turbulence in the postshock 
cooling layer is not available now. Thus we considered both 
limiting cases described above and show the uncertainty intro- 
duced by the lack of data concerning the large scale turbulence 
properties in Fig. 1. To avoid overestimating the nonthermal 
emission, we used a conservative minimum value of E c = 120 
keV for the radiative shock structure described above, although 
one could expect Ec ~ 20 keV in the most favorable case. 

Region IV. The density and temperature in the radiative shell 
depend on the preshock density, magnetic field and abundances 
(Shull & McKee 1979; Chevalier 1999). From the model 
of Chevalier (1999) for the above preshock parameters, we 
adopted n w 275 cm" 3 (N H wl.6x 10 20 cm" 2 ) and T ~ 10 2 
K in the dense radiative shell. The mass of the dense radiative 
shell is <~ 1000 M Q and its magnetic field is — 6 x 10" 5 G. The 
ionization in the dense radiative shell is supported at the level 
of about a percent by shock accelerated particles penetrating 
into the shell (see section 3.5 for details). The phase velocity 
of the high frequency collisionless Alfven waves which are re- 
sponsible for electron resonant scattering is about 50 km s" 1 ; it 
is determined by the number density of the ions. The Coulomb 
loss rate of the electrons is then dominated by the neutral com- 
ponent. 

We numerically calculated the electron distribution function 
in regions (I-IV). The results were then used in calculations of 
the nonthermal emission. 

3. NONTHERMAL EMISSION FROM IC 443 

The SNR IC 443 (G189.1+3.0) is a very good candidate for 
testing our model because the interaction of the SNR with a 
molecular cloud was established from radio and infrared line 
observations (e.g., DeNoyer 1979; Mufson et al. 1986; Burton 
et al. 1990; van Dishoeck et al. 1993; Richter et al. 1995; 
Claussen et al. 1997). The high energy fluxes are somewhat 
better determined than in other cases. Although we concentrate 
on IC 443, we briefly mention the case of W44, where observa- 
tions of OH masers and CO clumps again provide convincing 
evidence of molecular cloud interaction (Claussen et al 1997; 



Frail & Mitchell 1998). 

To model the nonthermal emission from the shell of IC 443 
we integrated the local emissivities over the radiative shock 
structure (regions / -IV). As a parcel of gas containing ac- 
celerated nonthermal electrons evolves through the radiative 
shock structure including the dense shell, the local emissivity 
from bremsstrahlung and inverse Compton was calculated us- 
ing standard theory (see Appendix A for details). The elec- 
tron - ion bremsstrahlung emission differential cross sections 
(Bethe-Heitler formulae) were taken from Akhiezer & Berestet- 
sky (1957) taking into account the Elwert factor, which is im- 
portant for modeling the keV emission. The electron - elec- 
tron bremsstrahlung contribution to the 7-ray emissivity was 
calculated using cross sections derived by Haug (1975). For 
the inverse Compton emissivity of IC 443, we used the same 
description of the background photon field as Gaisser et al. 
(1998). The hard emission production cross sections used in 
our calculations are similar to those used recently by other au- 
thors (Asvarov et al. 1990; Sturner et al. 1997; Gaisser et al. 
1998; Baring et al. 1999). We also estimated the contribution 
of pion decays to the 7-ray emission using a model for the pro- 
ton component and the cross sections from Dermer (1986). The 
resulting emission in the EGRET regime was roughly 7 times 
below that from electron bremsstrahlung, which is consistent 
with the results of Sturner et al. (1997). We do not discuss the 
pion component further. 




E ph (keV) 



FIG. 1. — Broadband uF v spectrum of the shell of IC 443 (distance 1.5 
kpc) calculated from a model of nonthermal electron production by a radia- 
tive shock with direct injection of electrons from the thermal pool. The shock 
velocity is 150 km s _1 , the interclump number density is 25 cm 4 , and the 
interclump magnetic field is 1.1 X 10~ 5 G. The electron diffusion coefficient 
(see equation [2]) is k im = 1.1 X 10 20 cm" 2 s _1 , a = 1.0, E(p,) = 1 MeV, and 
£(p**) = 20 GeV (equation [3]). The radiative shell parameters are described 
in § 2. The two solid curves labeled by numbers 1 and 2 correspond to two 
limiting values of Ec'. Ec = 120 keV (curve 1) and Ec = 2 GeV (curve 2). 
Interstellar absorption of the radio spectrum is treated as described in the text. 
The extended inverse Compton emission from the whole remnant is shown as 
a dashed line. The observational data points are from Esposito et al. (1996) 
for the EGRET source 2EG 10618 +2234 and Erickson & Mahoney (1985) for 
the IC 443 radio spectrum. Upper limits for MeV emission are from OSSE 
observations (Sturner et al. 1997). Upper limits for 7-rays > 300 GeV are 
from Whipple observations of IC 443 (Buckley et al. 1998). 

The synchrotron emission of the electrons calculated using 
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the scheme described by Ginzburg (1979) can be compared to 
the radio observations and to the upper limits at higher frequen- 
cies. We also included free-free absorption of low frequency 
radio waves (dotted line in Fig. 1). The synchrotron losses of a 
relativistic electron were included following Ginzburg (1979); 
bremsstrahlung energy losses are relatively unimportant in our 
case. 

3.1. Energetic Nonthermal Emission from the Shell 

The emission was integrated over the regions of the radia- 
tive shock in order to model the spectrum of the entire shell. 
Region IV is the dominant contributor to the synchrotron and 
bremsstrahlung emission. The results in Fig. 1 can be com- 
pared to the ASCA and EGRET observations of SNRs, as well 
as to the significant upper limits that have been set at TeV ener- 
gies by Whipple observations (Buckley et al. 1998) and from 
the preliminary analysis of CGRO OSSE observations presented 
by Sturner et al. (1997). We tried different parameter sets to fit 
the observations. 

In Fig. 1, vF v fits to the observations of IC 443 are presented 
for the model described in the previous section. The preshock 
density is taken to be n = 25 cm" 3 and the assumed distance to 
IC 443 is 1.5 kpc (e.g., Fesen & Kirshner 1980). Interstellar 
photoelectric absorption of the X-ray spectrum is accounted for 
with a line-of-sight column density Nh = 2 x 10 21 cm -2 using 
the cross sections from Morrison & McCammon (1983). 

For v S7 = 1 .5 in the diffusion model described by equation 
(2), we used k III0 = 1.1 x 10 20 cm 2 s -1 , a = 1.0, E(p t ) = 1 
MeV, and E(p^) = 20 GeV, which are compatible with the- 
oretical estimates of resonant wave-particle interactions. The 
corresponding value r\ k> 30. While in principle rj can be deter- 
mined from observations of heliospheric shocks, single space- 
craft observations do not provide direct values and are subject 
to ambiguous interpretation. Values of rj ~ 30 are consistent 
with modeling of the highly oblique solar wind termination 
shock (Ellison et al. 1999), but smaller values (rj ~ 3- 10) 
have been inferred for interplanetary traveling shocks (Baring 
et al. 1997). On the other hand, very weak scattering (rj <~ 100) 
has been inferred from observations near corotating interaction 
regions (e.g., Fisk, Schwadron, & Gloeckler 1997) and in the 
solar wind when the interplanetary magnetic field is nearly ra- 
dial (e.g., Mobius et al. 1998). We use rj w 30 as a conservative 
value. If strong MHD turbulence is present in the upstream 
shock region, the rj values could be closer to 1. This might 
allow an electron acceleration model similar to that described 
earlier, but for slower shock velocities ( < 100 km s" 1 ). 

The diffusion propagation model is not unique. We obtained 
a similar fit to that presented in Fig. 1 for a diffusion model 
with k im = 10 20 cm 2 s-\ a = 0.7, E(p*) = 1 MeV, %,) = 
20 GeV. As mentioned in §2, there is an important uncertainty 
concerning the lack of a quantitative model for the generation 
of large scale MHD turbulence in the postshock cooling layer. 
The uncertainty affects the hard X-ray spectral calculations. In 
view of this, we calculated the expected bremsstrahlung emis- 
sion for two limiting cases covering the range of uncertainties. 
Curve 1 in Fig. 1 corresponds to fully developed Alfvenic large 
scale MHD turbulence in the postshock cooling layer (E c =120 
keV). Curve 2 corresponds to the case of a lack of large scale 
MHD turbulence (E c = 2 GeV). While the GeV bremsstrahlung 
emission (as well as synchrotron and inverse Compton emis- 
sion) is similar in the 2 cases, the hard X-ray emission is sen- 
sitive to the uncertainty. We also obtained a similar fit to that 
presented in Fig. 1 for a model with v S7 = 1.1 and a diffusion 



model described by k um = 5 x 10 19 cm 2 s" 1 , a = 1.0, E(p*) = 1 
MeV, and E(p**) = 10 GeV. The model with « 1.1 requires 
a somewhat higher power conversion efficiency than that with 
v S7 w 1.5 (see §3.6), but this value is preferred for the IC 443 
radiative shock velocity (Chevalier 1999). 

The model of electron injection from the cloud thermal pool 
described above is based on the ionization structure of a ra- 
diative shock, which is very sensitive to the shock velocity if 
v S7 ~ 1 (Shull & McKee 1979; Hollenbach & McKee 1989). 
In this respect, it is instructive to consider a model in which 
high energy electrons are not related to freshly injected parti- 
cles, but are reaccelerated cosmic ray electrons (Blandford & 
Cowie 1982; Chevalier 1999; Cox et al. 1999). To model 
that case we suppose that the far upstream electron flux is just 
that observed for Galactic cosmic ray electrons near Earth, but 
extrapolated back to the energy E crm with the same slope as ob- 
served above a GeV. The calculations of vF v for this case are 
shown in Fig. 2. We consider the radiative shock structure with 
the parameters described in §2 (the same as were used for Fig. 
1). The Galactic cosmic ray electron flux in the far upstream 
region was taken to be 3x lO^^/GeV)" 3 e* s" 1 cm" 2 GeV -1 
sr" 1 (e.g., Berezinsky et al. 1990) for E > E crm . We assumed 
a flat cosmic ray electron spectrum below E crm , corresponding 
to flattening due to Coulomb losses in the interstellar medium. 
The maximum energy of cosmic ray electrons reaccelerated by 
the radiative shock is E m ^ 10 GeV (see equation [3]), as in the 
case shown in Fig. 1. The energy E crm was considered as a free 
parameter here. 

10 3 f , , , , , , a 
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FIG. 2. — Broadband vF v spectrum of the shell of IC 443 calculated for 
the model of Galactic cosmic ray electron reacceleration by a radiative shock. 
The basic shock parameters are as in Fig. 1 . The low energy cut-off of the 
preexisting Galactic cosmic ray electron spectrum in the cloud is 3 MeV for 
curve 1, 8 MeV for curve 2, and 50 MeV for curve 3. The inverse Compton 
emission from the whole remnant is shown as dashed lines. 

We do not distinguish between electrons and positrons be- 
cause there is no difference between them in the MHD shock 
acceleration process considered above and most of the emis- 
sion from high energy leptons involved in our model remains 
the same for both kind of particles. Because of interactions with 
the radiative shock structure, accelerated cosmic ray electrons 
could provide a good fit for both synchrotron radio and EGRET 
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emission if E crm < 10 MeV for the radiative shell parameters 
described above, but fall short for E crm « 100 MeV. 

The relatively low energies, E crm ^10 MeV, required for the 
IC443 shell are a potential problem for this model. The radio 
synchrotron radiation from the Galaxy indicates that the cos- 
mic ray electron spectrum flattens to N(E) oc E~ 2 below en- 
ergies <~ 1 GeV (Webber 1983). If we use such a spectrum 
in our model, the nonthermal emission is significantly below 
that observed. However, there are uncertainties about betatron- 
type acceleration in the shock compression, about the magnetic 
field strength, and about the volume of the emitting region. 
Duin & van der Laan (1975) did successfully model the IC 443 
radio flux with a cosmic ray compression model. We found 
that models with a higher magnetic field and/or larger volume 
could reproduce the radio flux, but that the 7-ray fluxes ob- 
served by EGRET could not be achieved. Cox et al. (1999) 
have recently concluded that cosmic ray compression by radia- 
tive shocks could provide a good fit to radio and gamma-ray 
emission from W44 if the cosmic ray electron population in 
the vicinity of the SNR is higher than that of measured in so- 
lar neighborhood. They used a somewhat higher magnetic field 
(B ~ 200 ^G) in the radiative shell of W44 than the 60 p,G field 
we assumed for the IC 443 shell. 

The acceleration time is a constraining parameter in this 
model because of the relatively low shock speed v S7 <~ 1 .5 in 
the dense medium. The test particle shock acceleration time 
t a (E) can be estimated from the equation 

u= ^_r(m^y± ( 5 ) 

vi-vnJ po V v/ v„ ) p 

where the velocities v, are measured in the shock rest frame 
(e.g., Axford 1981). For the radiative shock parameters de- 
scribed above, the acceleration time for E > 10 GeV is about 
3,000 years, for 77 ^ 30 (note that t a oc 77), consistent with the 
age of IC 443 estimated from the hydrodynamical model by 
Chevalier (1999). On the basis of Ginga observations of IC 
443, Wang et al. (1992) suggested an age of about one thou- 
sand years and an association with the supernova of AD 837. 
The reasons for the low age are related to a thermal interpre- 
tation of the observed X-ray emission above 10 keV. We shall 
argue below for the nonthermal origin of the hard emission in 
a radiative shock model. A distinctive feature of the radiative 
shock SNR shell model with direct injection of the electrons 
from the thermal pool is a hard spectrum of nonthermal emis- 
sion extending from keV to GeV energies (Fig. 1). There is 
also a hard emission component due to inverse Compton emis- 
sion of relativistic electrons having a scale comparable to the 
scale of the whole radio image of the remnant. 

3.2. Radio Emission from the Shell 

Observations of the radio emission from IC 443 were per- 
formed with different instruments over the last 35 years (see, 
e.g., Erickson & Mahoney 1985; Green 1986; Claussen et al. 
1997 and references therein). The 151 and 1419 MHz radio 
maps (Green 1986) reveal that the most intense flux comes from 
the shell of the remnant. We suggest that the eastern part of 
the shell is a radiative shock which has a high magnetic field 
and electron flux, dominating the synchrotron radiation from 
the remnant. 

We calculated the synchrotron radiation from the remnant us- 
ing the local emissivity given by equation (A17). The integrated 
spectrum from the remnant is dominated by the emission from 



the radiative shell (having magnetic field B~6x 10~ 5 G) if 
the magnetic field in the central parts of the remnant (region 
V) is moderate, B < 4x 10~ 6 G. For modeling the shell emis- 
sion, we took into account the interstellar free-free absorption at 
low radio frequencies, assuming a column density to IC 443 of 
Nh ~ 2 x 10 21 cm" 2 and an average interstellar medium temper- 
ature T ~ 10 3 K. This column density is similar to that adopted 
for the recent analysis of ASCA data by Keohane et al. (1997) 
and is about four times less than the value used earlier by Petre 
et al. (1988) and Wang et al. (1992). We assumed that the 
ionization fraction is ~ 15%. The free-free absorption inside 
the cool radiative shell is not very important above 50 MHz be- 
cause the emission measure is about 1 cm" 6 pc. The solid and 
dot-dashed lines on Fig. 1 show the total synchrotron radia- 
tion from IC 443 (dominated by shell emission in our model) 
expected at the Earth compared to the observational data (Er- 
ickson & Mahoney 1985). Both of our models are illustrated: 
direct injection of thermal electrons (Fig. 1) and reaccelera- 
tion of preexisting cosmic ray electrons by a radiative shock 
(Fig. 2). A radiative shock with total compression ratio 8, = 
4.4, including a shock precursor and subshock (S s = 3.3), was 
considered. An important feature of the radio observations of 
IC 443 and W44 is that their spectral indices are relatively flat, 
implying electron energy spectra flatter than those produced by 
standard strong shock acceleration in the test particle limit. In 
our model, the calculated radio spectrum fits that observed from 
the IC 443 shell because of the large shock compression ratio 
5, s=s 4.4 and to a lesser extent because of the effect of second or- 
der Fermi acceleration discussed recently by Ostrowski (1999). 

3.3. Nonthermal Emission from a Shocked Clump 

Molecular clouds show complex structure with an ensemble 
of dense massive clumps embedded in the interclump medium 
(Blitz 1993). We considered the possibility of radiative shock 
interaction with dense molecular clumps, resulting in high en- 
ergy emission. In IC 443 and W44, those regions do not appear 
to be significant sites of GHz radio continuum emission, but 
the clumps may contain most of the cloud mass and are possi- 
ble sites of hard X-rays and MeV 7-ray emission. The inter- 
action of a radiative shell with a molecular clump generates a 
dense slab bounded by two shocks (Chevalier 1999). The for- 
ward shock propagating into the clump is strong (see also the 
simulations of Jones & Kang 1993). Since the clump size is 
typically smaller or comparable to that of the radiative shock 
layer, it is embedded in a pool of nonthermal particles cre- 
ated by the radiative shock according to the scenario described 
above. Accelerated particles as well as ionizing radiation from 
the hot interior regions of the SNR should provide ionization 
of the external portions of a dense clump at the level > 1%. 
This is sufficient for a shock propagating into the dense clump 
with n <~ 10 4 cm" 3 and B <~ 0.1 mG (e.g., Claussen et al. 
1997) to be collisionless. The percent ionization level provides 
Alfven Mach numbers M a > 1 and weak damping of collision- 
less Alfven waves up to wavelengths ~ 10 10 cm for a flow with 
velocity > 25 km s" 1 in the clump. Then the injection and 
acceleration processes considered above would be efficient and 
the ionization in the preshock region would be maintained self- 
consistently at least at a percent level even for a shock propa- 
gating inside the dense clump. 

The gas ionization by UV radiation in the radiative shock 
precursor is sensitive to the shock speed if v sl < 1.1 (Shull & 
McKee 1979; Hollenbach & McKee 1989). This implies that a 
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shock having v$i < 1 propagating in a dense clump would be 
able to accelerate electrons only below E m < 1 GeV (see equa- 
tion [3]). We shall estimate below also the gas ionization rate 
due to accelerated nonthermal particles. The lifetime of a pre- 
existing electron of Me V energy against ionization losses in the 
interclump medium (n ~ 10 cm -3 ) is about 10 4 years. Unless a 
nearby source of fresh accelerated electrons is present, Galac- 
tic cosmic ray MeV electrons are unlikely to be present in the 
cloud. Direct injection of electrons from the thermal pool of the 
weakly ionized precursor would be a dominant source of shock 
accelerated electrons up to the cut-off energy E m . Because of 
a lack of accelerated GeV electrons in that case (v S7 < 1) the 
source would not appear as a bright radio and hard EGRET 7- 
ray emission feature. Nevertheless, these shocks might mani- 
fest themselves as hard X-ray and MeV 7-ray emission sources 
with vF v oc v a (a > 0.5) up to a few MeV and with a soft spec- 
trum at higher energies. A comparison of Fig. 3 with Fig. 1 
shows that the hard X-ray flux at 7 keV and above could be 
dominated by the emission from the shocked clumps, which 
is consistent with ASCA observations of IC 443 presented by 
Keohane et al. (1997). The hard X-ray nonthermal emission 
flux may be variable on a timescale of years. A relatively high 
MeV emissivity is a distinctive feature of shocked clumps, but 
they are difficult to observe. They must be nearby (within a 
few hundreds of parsecs) to be detected with COMPTEL/OSSE 
CGRO. 

For a forward shock of velocity v S7 < 1 propagating into a 
dense clump of number density n <~ 10 3 cm" 3 with upstream 
ionization level of a few percent (at <~ ki[E m ]/vs distance), the 
upper cut-off energy E m < 0.5 GeV. This is about 40 times less 
than that for a radiative shock in the interclump medium. The 
spatial diffusion k® below scales as w,A. Since the scale of 
a postshock radiative cooling layer oc n" 1 and A is a fraction of 
that, we have kmo ~ 3x 10 18 cm 2 s" 1 in the clump case. 

The effective compression of the nonthermal electrons as 
well as the magnetic field by the system of two shocks bounding 
the slab could be high. The magnetic field in the clump is ex- 
pected to be high. Claussen et al. (1997) derived line-of-sight 
magnetic field strengths of ^0.2 mG, remarkably uniform on 
the scale of several parsecs. Relativistic electrons accelerated 
by a shock propagating into a clump should produce radio emis- 
sion with a flat spectrum up to some hundreds of MHz, which 
is lower than the maximum radio frequency from the shell be- 
cause of the lower E m in the clumps. The emission measure of a 
clump interacting with a radiative shock could be EM < 1,000 
cm" 6 pc, which is much higher than that of the shell. Thus, the 
internal free-free absorption could be important here providing 
a flatter radio spectrum in the 100s of MHz regime. This is pos- 
sibly a reason for the observed spatial variations of the radio 
spectral index in IC 443 (e.g., Green 1986). 

We present in Fig. 3 the calculated vF v spectrum of nonther- 
mal emission from a radiative shock - molecular clump interac- 
tion region. The forward shock velocity was high - 100 km s" 1 , 
the initial number density in the clump 1,000 cm" 3 , and the 
shocked clump radius 0.5 pc. The shock compression ratio was 
the same as that in the radiative shell model: S s w 3 and S, w 
4.4. The ionization of gas in the shock upstream was about 5% 
(at the distances ~ ki[E m ]/v s ). To account for some uncertainty 
in the gas ionization, we present the radio spectra for two pos- 
sible values of the emission measure EM of the shocked clump: 
EM = 1 10 cm" 6 pc (solid line) and EM = 900 cm" 6 pc (dashed 
line). The diffusion coefficient was knio ~ 3x 10 18 cm 2 s" 1 , 



E(p^) = 1 MeV. We illustrate the results for two possible max- 
imal energies of the electrons (depending on the time of shock 
- clump interaction). Curve 1 in Fig. 3 corresponds to E m = 
0.5 GeV and E m = 0.05 GeV for curve 2. The synchrotron ra- 
dio spectrum for curve 2 is below the scale of the plot. The 
radio spectrum corresponding to curve 1 is very hard because 
of internal absorption in the shocked clump medium. 
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FIG. 3. — Broadband vF v spectrum for a model of nonthermal electron 
production by a radiative shock interacting with a molecular clump (distance 
1.5 kpc). The model spectra are calculated for the case of a shock velocity 
in the clump of 100 km s" 1 , a number density of 10 3 cnr 3 , and a magnetic 
field strength of 2 X 10~ 4 G. The electron diffusion coefficient (see equation 
[2]) k Im = 3 X 10 18 cm" 2 s" 1 , a = 1.0, and £(/>*) = 1 MeV. The maximum 
energy of accelerated electrons is E m = 0.5 GeV (curve 1) and E m =0.05 GeV 
(curve 2). The model spectrum for the case of a 30 km s" 1 shock in a clump 
with number density 10 4 cnT 3 is shown by curve 3 (see the text). Interstellar 
absorption of the radio emission is as in Fig. 1 . The internal absorption at the 
shocked clump is shown for EM f» 110 cm pc (solid line radio spectrum) 
and EM f» 900 cnr 6 pc (dashed line radio spectrum). The dot-dashed curves 
are the shell emission spectra from Fig. 1 given for comparison. 

We presented above the results for nonthermal emission of 
a strong radiative shock of v S7 <~ 1 propagating into a molec- 
ular clump. For the case of radiative shock interaction with 
the molecular clumps in IC 443, lower shock velocities, V57 ~ 
0.3, are more realistic (e.g., Cesarsky et al. 1999; Chevalier 
1999). The preshock ionization may be dominated mostly by 
nonthermal particles in such a shock because UV radiation is 
inefficient. A self-consistent model of MHD collisionless shock 
propagation into a dense clump requires simultaneous simula- 
tions of the MHD turbulence spectral properties and the non- 
thermal particle kinetics. Instead, we suppose that a collision- 
less shock transition exists for vsrj ~ 0.3 and n = 10 4 cm" 3 and 
assume the same diffusion model as described in §2. We calcu- 
lated the upstream gas ionization, finding that the ionized frac- 
tion x > 0.9 may hold up to depths about 10 15 cm" 2 due to 
shock accelerated electrons. The Mach number of the shock 
with V57 <~ 0.3 is typically below 5 because of heating of up- 
stream gas up to 10 4 K by accelerated particles and the com- 
pression ratio is lower than 4. The maximum energy of ac- 
celerated electrons in that case is typically below an MeV. We 
illustrate the calculated vF v spectrum of nonthermal emission 
from a 30 km s" 1 velocity shock (5 t = 3 ) interacting with a 
molecular clump of density 10 4 cm" 3 by curve 3 in Fig. 3; this 
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model corresponds to kjjjo <~ 7x 10 16 cm 2 s" 1 and £(/?*) = 10 
keV. Such collisionless shocks, if they exist, would provide a 
softer spectrum of nonthermal radiation than that for high ve- 
locity shocks with v sl > 1. A more comprehensive study is 
required to model the nonthermal emission from the low veloc- 
ity shocks in detail. A search for hard X-ray emission correlated 
with molecular emission predicted by the simplified model de- 
scribed above is possible with high resolution instruments like 
Chandra and XMM. 

3 .4. Hard X-ray Emission 

IC 443 was a target of X-ray observations with HEAO 1 (Pe- 
tre etal. 1988), Ginga (Wang et al. 1992), ROSAT (Asaoka 
& Aschenbach 1994) and ASCA (Keohane et al. 1997). ASCA 
GIS observations discovered the localized character of the hard 
X-ray emission (Keohane et al. 1997). Most of the 2-10 
keV GIS photons came from an isolated emitting feature and 
from the southeast elongated ridge of hard emission. The ridge 
is also coincident with the 95% confidence error circle of the 
EGRET CGRO source (Esposito et al. 1996), while the iso- 
lated feature is outside the circle (Keohane et al. 1997). The 
integrated 2-10 keV flux observed by ASCA GIS was (5±l)x 
10~ n erg cm" 2 s" 1 , representing more than 90% of the total flux 
from the SNR (Keohane et al. 1997). The hard X-ray ridge, 
as defined by Keohane et al. (1997), is only a part of the more 
extended radio shell discussed in the previous section. The flux 
density at 7 keV from the isolated emitting feature was about 
4x 10~ 5 photons cm -2 s" 1 keV" 1 and that from the ridge was 
2x 10~ 5 photons cm" 2 s" 1 keV" 1 . HEAO A-2 measured a flux 
at 2-10 keV of (7±l)x 10" 11 erg cm" 2 s" 1 (Petre et al. 1988), 
while Ginga measured a 2-20 keV flux of about 9x 10" 11 erg 
cm" 2 s" 1 (Wang et al. 1992). These numbers are consistent 
with the assumption that the integrated flux is dominated by the 
extended soft component coming from the central parts of the 
remnant. 

The soft X-ray 0.2-3. 1 keV surface brightness map of IC 443 
from the Einstein Observatory (Petre etal. 1988) shows bright 
features in the northeastern part of the remnant. The presence 
of nearly uniform X-ray emission from the central part of the 
remnant is a characteristic feature of mixed-morphology SNRs 
(Rho & Petre 1998). It might be due to the effect of thermal 
conduction (Cox et al. 1999; Shelton et al. 1999), although 
Harms et al. (1997) argued for an alternative radiative shock 
model for W44, another mixed-morphology SNR. That feature 
in our model corresponds to the emission from hot (T < 10 7 K), 
low density gas interior to the shock (region V). 

In the radiative shock model with direct injection from the 
thermal pool and second order Fermi acceleration in postshock 
cooling layer, most of the hard X-ray emission ( > 7 keV) 
comes from the shell - radiative shock structure with a rising 
vF v X-ray spectrum (curve 1 in Fig. 1 for the shell and curves 
1 and 2 in Fig. 3 for a shocked molecular clump). In addition to 
the shell, there is a more extended source of hard X-rays from 
inverse Compton emission. It is the dominant source of hard 
X-rays in the scenario with cosmic ray electron reacceleration 
(Fig. 2) and in the case of a lack of second order Fermi ac- 
celeration in the radiative shell (see curve 2 in Fig. 1). In the 
scenario with reacceleration of cosmic ray electrons, it is im- 
possible to obtain the 7 keV flux observed by ASCA from the 
IC 443 ridge. Thus, a resolved hard X-ray image of the IC 443 
remnant may be very informative in determining the electron 
acceleration scenario. 

The isolated emitting feature in the southern part of IC 443 



detected by ASCA and brightest around 7 keV, has a flat spec- 
trum low frequency radio continuum source (Green 1986) and 
is close to a region of excess H2 emission. High spatial resolu- 
tion observations of this isolated emitting feature are required 
to decide between a low-luminosity pulsar nebula and shock in- 
teraction models (e.g., radiative shock interaction with a dense 
molecular clump, Fig. 3). Observations of the hard emission 
with the sensitivity and resolution available with BeppoSAX, 
Chandra XRO, XMM and GLAST would be required to con- 
strain the extension and spectral properties of the hard X-ray 
emission from IC 443. 

3.5. Gas Ionization by Nonthermal Electrons 

Nonthermal particles accelerated by a MHD collisionless 
shock wave provide an efficient ionizing agent. We calculated 
above the nonthermal emission generated by accelerated ener- 
getic electrons; one can also calculate the ionization produced 
by the same electrons in the molecular cloud. This might con- 
nect the high energy observations with radio, IR and optical 
emission from shocked atomic and molecular gas. To estimate 
the gas ionization due to accelerated electrons in the radiative 
shock structure we used the electron spectra calculated with 
the kinetic model of electron acceleration and propagation de- 
scribed in §2 and electron ionization cross sections from Ap- 
pendix B. In Fig. 4, we present the calculated ionization rate ( e 
and estimated ionized fraction x in the shell of a radiative shock 
due to accelerated electrons for the case of vsi = 1.5 for both 
models of large-scale turbulence in the postshock cooling layer 
considered in §3.1. 
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FIG. 4. — Postshock ionization rate due to nonthermal electrons Q (lower 
panel) and the estimated ionization fraction x (upper panel) versus total H col- 
umn density. The solid lines correspond to the model with Eq = 120 keV (see 
curve 1 in Fig. 1) and the dotted lines to Ec = 2 GeV (see curve 2 in Fig. 1). 

There should also be contributions to x from accelerated nu- 
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cleons and from UV radiation which are not included in Fig. 
4. We accounted for only radiative recombination, assuming a 
small molecular fraction in the radiative shell at column densi- 
ties below 2 x 10 20 cm" 2 . 

The solid curves in both panels of Fig. 4 correspond to the 
case of fully developed Alfvenic turbulence in the postshock 
cooling layer, which provides efficient second order Fermi ac- 
celeration making Coulomb losses relatively unimportant above 
E c = 120 keV. This case corresponds to curve 1 in Fig. 1 for 
emission from the shell with substantial hard X-ray emission. 
The dotted lines in Fig. 4 correspond to the case of a lack of 
large scale MHD turbulence where the Coulomb losses over- 
come the second order Fermi acceleration up to E c = 2 GeV. 
The nonthermal emission expected in this case is illustrated by 
curve 2 in Fig. 1. The radiative shell ionization by nonthermal 
electrons is sensitive to the turbulent structure of the postshock 
cooling layer and there is a correlation between the hard X-ray 
spectrum and the ionization structure of the radiative shell. We 
shall consider the effect of nonthermal particles on the ioniza- 
tion structure of a shocked clump elsewhere. 

3.6. Energy in Nonthermal Electrons 

The energy of the electron component for the parameter set 
used to compute the spectra in Fig. 1 is substantial. The power 
required for electron acceleration and maintenance for curve 1 
in Fig. 1, ~ 1.5 x 10 37 erg s -1 , is high because of the strong 
Coulomb losses of keV electrons in the dense ionized medium, 
but it is consistent with the estimate of the total radiative losses 
from IC 443, dominated by infrared emission (Mufson et al. 
1986; Burton et al. 1990). A similar power is required for 
both models 1 and 2 in Fig. 3, where radiative shock - molec- 
ular clump interaction is illustrated. For curve 2, the power 
requirements are less, ~ 6 x 10 36 erg s" 1 . The implied effi- 
ciency of power conversion from the MHD shock flow (which 
can be estimated as 3 xlO 38 erg s" 1 ) to the nonthermal elec- 
trons is about 5% in these cases. The model with v$i ~ 1 would 
provide emission similar to that shown in Fig. 1, but requires 
about three times higher efficiency. This is higher than the elec- 
tron acceleration efficiency estimated from GeV regime cosmic 
rays observed near the Earth, but evolved SNRs interacting with 
dense molecular gas probably cannot be considered as the main 
source of the observed Galactic cosmic rays because of the rel- 
atively low maximum energies of accelerated particles. 

The pressure of the nonthermal electron component down- 
stream of the shock (region IV) is V e ~ 7x 10 -11 erg cm" 3 . 
The magnetic pressure in the dense shell is about 1.5 x 10~ 10 
erg cm -3 , which is higher than the thermal gas pressure. The 
uniform (~ 60 fiG) and stochastic magnetic field components 
dominate the total pressure in the shell. 

Let us consider the energy requirements for the scenario with 
electron injection from preexisting cosmic rays (Fig. 2). The 
nonthermal electron pressure in this case is V e ~ 5.2 x 10" 11 
erg cm -3 (for E crm < 8 MeV), similar to that for the model with 
injection from the thermal pool, because it is dominated by GeV 
particles. An advantage of this scenario is that much less power 
is required. It is about 4 x 10 35 erg s" 1 , which is less than 10% 
of that for the scenario with injection from the thermal pool. 
The difference is because Coulomb losses are unimportant for 
the high energy electrons involved in that scenario. However, a 
potential problem is the relatively low E crm required. The life- 
time of a 20 MeV electron against ionization losses in a neutral 
medium of number density <~ 1 cm" 3 is about 10 6 years. This 
implies that the source of MeV cosmic ray electrons should ex- 



ist within about 100 pc of the molecular cloud if the diffusion 
coefficient is ~ 3x 10 26 cm 2 s" 1 at these energies. 

4. DISCUSSION AND FUTURE PROSPECTS 

For individual supernova remnants, the ambient density is 
an important parameter that can be estimated from multiwave- 
length observations (Chevalier 1999). Another important pa- 
rameter for modeling the nonthermal particles in a SNR in 
a molecular cloud is the collisionless MHD turbulence spec- 
trum, particularly in the radiative shock cooling layer. A sub- 
stantial level of MHD collisionless turbulence could overcome 
Coulomb losses in the dense plasma downstream from the ra- 
diative shock. An accurate model of MHD turbulence in the 
postshock cooling layer is not available now. Thus we consid- 
ered both limiting cases described in the § 2 and show in Fig. 
1 the uncertainty introduced by the lack of data concerning the 
MHD turbulence properties. Coulomb losses are important in 
the postshock cooling region only for electrons with energies 
below Ec ~ 20 keV if the large scale turbulent velocity has a 
substantial longitudinal component of <~ 20 km s -1 . A higher 
value (E c > 8 MeV) is expected in the absence of a longitudi- 
nal component of large scale turbulence in the postshock cool- 
ing layer because the Coulomb losses are overcome by resonant 
interaction with Alfven waves. Since the Alfven velocity is rel- 
atively low (~ 5 km s" 1 ), the value of Ec is much higher in that 
case. We used a conservative minimum value of E c = 120 keV 
for the radiative shock structure described above, although one 
could expect Ec ~ 20 keV in the most favorable case. The non- 
thermal emission of SNR in a molecular cloud (especially in 
the hard X-ray regime) and ionization structure of the radiative 
shock are sensitive to the MHD turbulence model. Spatially 
resolved observations of the nonthermal emission from SNRs 
may be able to constrain models of MHD turbulence. 

All of the SNR candidates from the EGRET list: IC 443, 
7 Cyg, W44, and Monoceros (Esposito et al. 1996) are old 
remnants interacting with molecular clouds. The 7-ray emis- 
sion from electrons accelerated by the radiative shock calcu- 
lated for the IC 443 parameters presented in Figs. 1 and 2 is in 
good agreement with that observed by EGRET (Esposito et al. 
1996), as well as with the upper limits established by Whipple 
(Buckley et al. 1998). The high energy 7-rays ( > 50 MeV) 
from the radiative shell should be spatially correlated with the 
radio emission. The EGRET telescope detected an extended 
excess (95% confidence circle ~ 42' in radius) from the Mono- 
ceros SNR correlated with the 1 .42 GHz radio emission (Es- 
posito et al. 1996). This is in accordance with the model of 
a leptonic origin of GeV emission from extended SNRs inter- 
acting with clouds (de Jager & Mastichiadis 1997). GLAST 
(Gamma-ray Large Area Space Telescope) will be an excellent 
instrument for future 7-ray observations, because it will have 
the capability of spatially resolving the 7-ray emission from 
SNRs. 

Spatially resolved observations of W44, IC 443, 3C391 and 
some other mixed-morphology SNRs from the list given by Rho 
& Petre (1998) with BeppoSAX, Chandra XRO, XMM would 
be very valuable tools to test our model. The hard X-ray de- 
tector (HXD) aboard the forthcoming ASTRO-E mission could 
be used for observations of the hard 10-700 keV continuum 
with a field of view of 0.8° FWHM at 60 keV and 2.8° at 500 
keV from the extended hard X-ray emission from IC 443, W44 
and the Monoceros SNR predicted by the radiative shock model 
(see Fig. 1 and Fig. 3). Spatially resolved spectra are needed 
to distinguish the shell emission from the shocked clump emis- 
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sion. 

Molecular clumps interacting with moderately fast radiative 
shocks are also expected to be sources of hard X-rays and MeV 
7-rays (up to 100s of MeV). Nonthermal continuum radio emis- 
sion (100s of MHz) with a sharply rising vF v spectrum (Fig. 3) 
and a time dependent cut-off frequency is expected in the fast 
shock model. Due to potentially substantial internal free-free 
absorption, the spectrum of radio emission from a localized 
clump might constrain the ionized gas density in the clump. 
The MeV 7-ray spectrum of a localized clump may be resolved 
with the forthcoming 7-ray missions INTEGRAL (e.g., Schon- 
felder 1999; Winkler 1999) and GLAST. With an expected an- 
gular resolution about 12' FWHM, imager IBIS aboard INTE- 
GRAL would allow detection of hard X- ray ( > 50 keV) emis- 
sion from molecular clouds within ~ 1 kpc. A comparison 
of hard emission spectra with the radio spectrum can provide 
valuable information about the density and magnetic field in 
a clump. A possible variability on a timescale of a few years 
for clump hard emission and radio emission (on a longer time 
scale) could further constrain the model. 

For low velocity (below ~ 30 km s" 1 ) shocks interacting with 



a dense (~ 10 4 cm -3 ), magnetized (B < 0.5 mG) molecular 
clump, one may expect hard X-ray emission below the MeV 
regime correlated with the regions of molecular emission of 
shocked gas. Nonthermal radio continuum is not expected in 
the case of low velocity shocks because of a lack of accelerated 
relativistic electrons at sufficiently high energies. 

Finally, we note that the high density of energetic particles 
in the vicinity of a shock wave in a molecular cloud can affect 
the ionization and thermal properties of the gas. Rich molecu- 
lar spectra have been observed from IC 443 and 3C391 (Reach 
& Rho 1999). It will be interesting to see whether there is a 
signature of the presence of nonthermal particles that can be 
discerned from the molecular spectra. 
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APPENDIX 

APPENDIX A: RADIATIVE PROCESSES 

Given the electron distribution function, the high energy emission flux J{E^) at a distance R from the source can be calculated 
from the equation: 

1 f nr dn^(E 1 ^f) 



J ^=4^J dv ^r^ (A1) 

v 

where dn 1 (E 1 ^f)/dt is the emissivity, which here includes bremsstrahlung, synchrotron radiation and the inverse Compton effect. 
We generally assumed a pure hydrogen composition in our calculations. 

High energy particles penetrating through a partially ionized medium produce photons due to interactions with atoms, electrons 
and ions. We assumed a Maxwellian distribution for the ambient matter and in most cases took the target ambient particles to be at 
rest. Then, the emissivity can be calculated from 



dn~f(E 7 ) 
dt 



■ At: 



(00 00 ^ 

J dE e N p a e p(E 7 ,E e )J(E e )+ J dE e N e a ee (E j7 E e )J(E e ) 

\E-, £„„■„ / 



(A2) 



where J(E e ) is the electron flux measured in units s" 1 cm" 2 keV" 1 sr" 1 , cr ep and <j ee are differential cross sections of one photon emis- 
sion due to electron-proton and electron-electron interactions integrated over angles, and E min is the minimum energy of incoming 
electron that could generate a photon with energy E 1 . For e-p interactions, £,„,„ = £ 7 is a good approximation. The correction due to 
finite proton mass is less than the uncertainty in the cross section. In the case of e-e interactions, we used the following equation for 
F • 

'-'mm ■ 



Ery — E m l n ffl e C I 



2,tn e c -\-E m j n \f ' E m j n {E m i n + 2fw e c 2 ) 



The inverse bremsstrahlung contribution is negligible unless there is a very high (far in excess of 100) ratio of protons to electrons. 
We used the following cross section approximations in our calculations: 

e-e bremsstrahlung. 

Haug (1975) has obtained a general formula for the e-e bremsstrahlung cross section, but it is more convenient to use simplified 
approximations. We used relativistic and non-relativistic asymptotic expressions, which were matched in the intermediate energy 
region (see Baring et al. 1999). 

In the relativistic limit, we used approximate formulae given by Baier et al. (1967) with a correction factor from Baring et al. 
(1999):0 

cr ee (E ei E 1 ) = (ai + <T2)A(e 1 ,"fe), (A3) 

4 Equations (A4) and (A6) are somewhat different from that given by Baring et al. (1999). M.G. Baring (private communication) made available the corrections 
to their original formulae. 



where 7, = (E e + m e c 2 )/m e c 2 , e 7 = E 7 /m e c 2 , 



4r 2 a 

0-i(7e,e 7 ) = 



1 + 



3 7,-1 



7,-1 



In 2(7,-1) 



7,-1-e- 



02 : 



16(l-e 7 + e 2 7 )ln(^)-^ + ^-4+4e 7 -8e 7 
-2(1 - 2e 7 ) ln(l - 2e 7 ) - ^ + 1- - 2 + 4e 7 ) 

( 4 -^ + 4^)ln(27,)-24-4 



r„ = e 2 /m e c 2 is the classical electron radius, a = e 2 /hc, and 



e < - 

e 7 ^ 2' 
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(A4) 



(A5) 



A(e 7 ,7 e ) = l- 



10(7,-1)3 /e 



3 7, + l V7< 



(A6) 



According to Baring et al. (1999), these formulae agree with those from Haug (1975) with an accuracy ~ 10% for electrons of 
energy ;> 5 MeV . 

In the non-relativistic limit, we used equations from Garibyan (1952) and Fedyushin (1952) with correction coefficients from 
Baring et al. (1999): 



4e 7 



(A7) 



where the photon energy is in the range < e 7 < j(7, - 1) and F(x) (0 < x < 1) has the form: 



F{x) = B{ le ) 
+ 



' (•(-,. v) 

7.t 2 3x 



(2-xf 

12(2-jc)-£- 



The correction coefficients are: 



2-x (2^F 

m ^ ,1,2 n rl x 10*7«ft(2 + 7«fc) 
2 l+x z (7g-l) 



(A8) 



(A9) 



Baring et al. (1999) found less than a 10% difference between equation (A7) and the expressions of Haug (1975) for electron 
energies E e <; 500 keV. 

e-p bremsstrahlung. 

In the non-relativistic limit, we took the e-p bremsstrahlung cross section from Akhiezer & Berestetsky (1957): 

,2 7 2_2„2 4jT 2 (iCi 



16r z Q Z z amic 2 



l ln [ Pl +P2 



P 3pj (e 27r C'-l)(l-e- 2 ^)e 7 \P1-P2J' 

which is valid if C,\ <C 1; here Q = aZ/(3i. 

For the relativistic limit, we used Born approximation formulae from Akhiezer & Berestetsky (1957): 



(A10) 



J ep ■ 



o />i (3 



p\ pipi 



ipTF2 + pJ7 2 (e i £ 2 +p i p 2 )+ 2p^ [ m ^~ m—^ + ^r 



Pi 



(All) 



where L = In ( ^^^ ) = 21n ( P ' + ^" m ' ), m = In (f^[) = 21n (^±f ), and m = 21n (^±f ). In the velocity range Za{f3? - 

/V), we included the Elwert factor from Koch & Motz (1959): f E = ^-e! P \-2^my 

In equation (Al 1), c = 1 is assumed. In the intermediate energy range, we matched equation (A10) and equation (Al 1). 

Synchrotron radiation. 

We used standard formulae for synchrotron radiation in our calculations (e.g., Ginzburg 1979). The emissivity can be written as: 



F dn-,{E-,) 
1 ' dt 

F(x) 



m e c 
00 



■jd X dEN(E,x)Hsm 2 (x)F^, 



xfK 5/3 (r])dr], 



(A12) 



12 

where \ is the angle between the electron velocity and the magnetic field, N(E,x) is the electron distribution function over energy 

and angle, v c = 3 ^"" ( c x) ) . and K 5 / 3 (x) is the McDonald function. Integrating equation (12) over angle under the assumption 
of a chaotic magnetic field orientation, we obtain the equation 



dt 2nm e c J \i> c 



(A13) 



where N(E,~~t) is the electron distribution function as a function of energy and position and H±_ = ^/2/3H(T > ) is the averaged magnetic 
field. 

Inverse Compton emission 

High energy electrons colliding with photons result in high energy photon production (e.g., Jones 1968; Gaisser et al. 1998; 
Sturner et al. 1997; Baring et al. 1999). If the electron energy E e <C ss £ 6 eV (E 7 is the energy of the emitted photon), 
then the formulae in the Thompson limit are (e.g., Ginzburg 1979): 

= V3ira T m e c 2 J dE ph ^^N ph (E ph Jt)J e (m e c 2 ^ ' ^L^j , (A14) 
where oj = -r- 

(£) 

is the Thompson cross section, N p h(E p h) is the background photon energy spectrum, and J e (E) is the incoming 



dn^iEy) 
dt 



electron flux. 

If the electron energy E > & ' %; eV, we must use the Klein-Nishina formulae (integrated over angles): 

m-y £ 7 (cV) 



' ^4?r / dE e dE ph N ph (E ph )a KN (E^,E e ,E ph )J(E e ), (A15) 



dt 

2lTr 2 - ..in _„>r: 



&KN(E 7 ,E e ,E p h) - 

ZphTe 



2q\xv{q)+\+q-2q 2 +^^ 1 



2(1+ 9 D 



(A16) 



where T = 4e ph ~/ e , q = 2—, 0<q<l. 

(7 e -e 7 )T 



In the above, e 7 = E y /m e c 2 , e p h = E p h/m e c 2 , j e is the electron Lorentz factor, and r„ = e 2 /m e c 2 is the classical electron radius. The 
interstellar background photon spectrum, N p h(E p h,R), was adopted from Mezger et al. (1982) and Mathis et al. (1983) taking into 
account infrared data from Saken et al. (1992) for the case of IC443 (see also Gaisser et al. 1998). 

Free-free absorption 

Radio waves propagating in the partially ionized thermal plasma are subject to thermal free-free absorption. For hv <C kT (the 
Rayleigh - Jeans regime), the free-free absorption coefficient oJJ (cm -1 ) can be calculated from 

ff 4e 6 / 2tt \ 1/2 , , 

Using an approximate formula for the Gaunt factor gff(y, T) from Rybicki & Lightman (1979) appropriate to our parameter range, 
we obtained a simplified expression for the optical depth t{v, T) 

rw 0.022 uiftT^EM, (A18) 

where the radio wave frequency ^100 is measured in units of 100 MHz, the plasma temperature in units of 100 K is 7^ 00, and the 
emission measure EM is in units of cm" 6 pc. 

APPENDIX B: ELECTRON IONIZATION 

The primary ionization rate ( e due to electron impact was calculated from the equation: 



00 

( e =4n J dEa eion (E)Je(E), (Bl) 

Emta 

where J e (E) is the electron flux measured in units s -1 cm -2 keV -1 sr -1 , a e i on is the differential cross section for hydrogen atom 
primary ionization by electron-atom interactions integrated over atomic electron states, and E m j„ is the minimum energy of the 
ionizing incoming electron. 
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In the regime below 1 keV, we obtained the following fit to the ionization cross section <7 e;o „ using experimental data from Fite & 
Brackmann (1958): 

aeion(E) = nalAi sxp(-bA 4 -b 3 A 5 -b 5 A 6 )/E(eVf 3 , (B2) 

where b = Ai/EieV) and E(eV) is the incoming electron energy measured in eV. The fitting coefficients are Ai=76.1 1, A2=14.34, 
A 3 =0.89, A 4 =3.82, A 5 =-2.55, andA 6 =4.5. a Q = h 2 /m e e 2 is the Bohr radius. 

At high energies ( > 1 keV), the Born approximation provides good accuracy (e.g., Mott & Massey 1965): 



eion — 



I 0K (K)dKdK, 



(B3) 



where 



k K (K) 



2tt2 10 k 



\k 2 K 



x exp 



K 2 + (fj 2 +K 2 )/3 



t i 4 + 2^ 2 (K 2 + n 2 ) + (K 2 -K 2 Y 

2fJ,K 



2a 

arctg. 



1 -exp(-27r/z/V 



(B4) 



u = Z/a , k(E) is the incoming electron wave number, k K is the outgoing electron wave number, and n is the emitted atomic electron 
wave number. In addition, K = \ k - k K \, K max = k+k K , and K min = k-k K . 

We found a convenient approximation providing a reasonably accurate fit (better than 10% above 1 keV and better than 5% above 
3 keV) to equations (B3) and (B4): 

a eio „(E) = ira 2 P l E(eVy P2 , (B5) 

where Pi=97.72 and P 2 =0.93. 

The gas ionization fraction x = rij/ (rij + n n ) (we consider here a hydrogen plasma) can be expressed through the ionization rate £ 
and recombination rate a r (T) as: 

x= (s 2 /4 + s) -s/2 (B6) 

where s = (/(a r n). We used the radiative recombination rates given by Spitzer (1978) and dissociative recombination rates from 
Schneider et al. (1994). 
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